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SUMMARY 


"—'An  aerothermochemlcal  analysis  of  chemically  active 
surfaces  immersed  In  a  high  speed  stream  of  an  oxidizing  gas  mixture 
is  presented.  The  analysis  applies  lu  turbulent  flow  over  blunt 
bodies  whose  surface  material  passes  into  the  gaseous  phase  without 
going  through  the  liquid  stste.  Thi*  report  -is  a  -  sequel— to  the 
previous  investigation1-1)  of  laminar  flow  under  the  same  environ¬ 
mental  condition*.  The  theory  is  developed  In  general  for  a  whole 
class  of  sublimators  and  is  then- applied  to  the  particular  case  of 
graphite  combustion  about  which  a  great  deal  of  fundamental  data  13 
available . 
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SECTION  1 


INTRODUCTION 


Thu  study  oi'  combustion  processes  requires  a  blending  of 
complex  chemical  and  tiuid  dynamical  phenomena.  Under  certain 
conditions  the  interactions  between  these  fields  can  be  minimized. 

For  example,  if  the  speed  of  chemical  reactions  is  rate  controlling, 
attention  can  be  focused  on  chemical  phenomena  with  little  emphasis 
on  the  accompaning  flow  conditions.  Here  the  most  difficult  problem 
to  be  faced  is  that  of  adequately  describing  the  chemical  kinetics 
Generally,  the  detailed  reaction  steps  are  r.ot  well  known,  and  even 
If  they  are  it  is  unusual  If  the  corresponding  Specific  reaction 
rate  coefficients  are  known  within  several  orders  of  magnitude. 

This  difficulty  is  often  circumvented  by  assuming  either  very  fast 
or  very  slow  reactions.  On  the  other  hand,  if  diffusion  and  convec¬ 
tion  are  rate  controlling,  prediction  of  combustion  phenomena 
depends  largely  upon  the  state  of  knowledge  for  the  particular 
viscous  flow  phenomenon  involved. 

Viscous  flows  can  be  classified  as  laminar,  turbulent  or 
trsnsitiun;  forced  convection  or  free  convection;  boundary  layer  or 
wake  flows.  Of  all  the  various  combinations,  laminar,  forced  con¬ 
vection,  boundary  layer  flow  lends  itself  most  easily  to  mathematical 
description.  Although  turbulent  flow  is  net  nearly  as  well  under¬ 
stood  in  detail,  turbulent  skin  friction,  mass  transfer,  and  heat 
transfer,  under  conditions  of  torced  convection  boundary  layer  flow, 
can  be  predicted  with  some  degree  of  success  with  the  aid  of  a  few 
empirical  constants  obtained  from  experiments.  It  certainly  does 
not  rank  last  among  the  various  possible  combinations  of  the  above 


I 


Ail  '•UBUCATICK  KO.  l'-'66 


INTRODUCTiTN 


phcntcuna .  However,  because  the  funoamentalr  rema.'.i:  rather  obscure 
any  attempt  to  make  predictions  Is  hold  In  rather  lew  esteoui  In  some 
circles.  Unfortunately , the  engineer  is  faced  with  turbulent  flow 
far  more  often  than  he  can  depend  upon  laminar  flow.  In  the  case 
under  consideration,  the  combustion  of  chemically  active  sublimating 
surfaces,  this  Is  even  more  likely  tc  be  true.  Both  the  high  sur¬ 
face  is:  perature  and  th*  '’■rfarp  notes  injection  should  tend  to 
destabilize  the  laminar  boundary  layer. 

The  laminar  case  was  investigated  by  the  author  and 
Ur.  D.  A.  Dooley,  also  of  Aeronutronic  Systems,  Inc.d)*  The  pre¬ 
sent  report  should  be  regarded  as  a  sequel  to  the  previous  work. 

Much  the  same  approach  la  uned  in  tiie  present  analysis.  In  the  gas 
phase,  diffusion  is  assumed  to  be  rate  controlling  so  that  the  chem¬ 
ical  reactions  chosen  to  represent  the  combustion  process  can  be 
considered  to  he  nearly  in  equilibrium.  At  the  surface  a  somewhat 
more  general  treatment  la  used.  The  solid  surface  combustion  and 
sublimation  reactions  are  considered  to  proceed  in  such  a  way  that 
the  products  are  swept  away  by  diffusion  and  convection  faster  then 
they  are  produced  when  the  flux  of  species  is  high;  while  if  the 
flu*  rate  is  low,  the  hackscattering  of  products  is  sufficient  to 
approximate  equilibrium.  This  amounts  to  allowing  either  chemical 
reaction  or  diffusion  to  be  rate  controlling  at  the  surface  under 
appropriate  circumstances. 

Ir.  its  boundary  layer  aspects  the  present  semi-empirical 
theory  uses  an  approach  similar  to  that  of  Van  Driest^)  O)  who 
allowed  for  variations  in  density  through  the  J>oundary  layer  in  such 
a  manner  that  the  results  reduce  to  the  von  Karman-Schoenherr 
relationship  for  incompressible  skin  friction.  It  has  been  pointed 
out  that  there  are  several  iurmulatlons  for  including  variable 
density  which  can  be  made  to  reduce  to  the  incompressible  results^) 
but  Van  Driest 's  physical  reasoning  seems  to  have  been  sufficiently 
good  so  that  his  results  agree  well  with  experiments.^  It  is 
hoped  that  an  extension  of  this  method  to  thecase  of  combustion  will 
meet  with  some  measure  of  success.  However,  until  experimental 
evidence  is  available  the  results  obtained  should  be  used  with 
caution. 


*  Since  publication  of  that  report  a  great  deal  of  interest  has 
been  stimulated  In  this  field.  Papers  by  Lees^)  and  by  Cohen, 
Bromberg  and  Llpkl s  (-'J  treat  similar  chemical  systems  with 
similar  results. 


2 


SRC  1 1DS  2 


T', !  i.N,  .iiKT.UAi;',  R^’.'TIOMS 


Although  prediction  ot  turbulent  boundary  layer  phenomena, 
within  the  present  state  of  knowledge,  must  be  Semi-ercpi rical  in 
nature,  analysis  of  the  fundamental  equations  helps  form  a  logical 
framework  for  physical  reasoning.  The  development  of  the  funda¬ 
mental  equations  resfc  on  the  conventional  assumption  that  the 
conservation  equations  to.  a  continuum  hold  instantaneously.  The 
instantaneous  value  of  each  of  the  unknowns  is  separated  into  the 
sum  of  a  mean  and  a  fluctuating  part.  Introduction  of  such  expres¬ 
sions  into  the  conservation  equations  for  reacting  gas  mixtures  and 
averaging  according  to  the  Reynolds  rules  for  tine  averagrs  yield 
average  conservation  equations  which  include  correlations  between 
fluctuating  components.  The  full  equations  which  result  are 
presented  in  Appendix  A  where  the  s impl i t icat ions  which  can  be  made 
prior  to  introducing  the  boundary  layer  assumptions  are  discussed. 

As  in  the  laminar  analysis,  it  is  assumed  that  diffusion  due  to 
temperature  and  pressure  gradients  is  negligible  compared  to  that 
due  to  concentration  gradients  and  that  all  diffusion  coefficients 
are  equal.  When  correlations  involving  fluctuations  in  molecular 
transport  terms  are  neglected  according  to  the  reasoning  of 
Appendix  A  and  boundary  layer  assumptions  are  made,  i.e., 

Vnsjt  and  K  g  small,  the  following  boundary  layer 
equations  result  in  the  mean  steady  state  for  blunt  body  coordinate 
systems : 
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In  this  and  the  following  equations  the  bars  denote  averages  and  the 
primes  denote  fluctuating  components.  The  quantity  K:  is  the  mass 
fraction  of  the  ith  molecular  species  while  (Af.  is  its  rate  of 
production  per  unit  volume. 

Overall  Continuity 

\(?7Fa,w) 

+  ^  ”  °  <2> 


where  &  is  unity  for  bodies  of  revolution  and  zero  for  two 
dimensional  bodies. 


The  atomic  mass  fractions,  »  are  related  to  the  molecular  mass 

fractions  through  the  definition: 


<A) 
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where  /2  ^ ^  represents  the  mass  of  atf.-t.iic  species  fa  per  unit 
n.ess  a!  molecular  species  i 

Vomcntum 


v  >-/ 


♦ 


S  ia/  . 


-  yn  u  ^ 


(  p  *(?*-)' /v‘  ) 

■ 


(6) 


where  ft.  and  are  the  density  and  velocity  at  the  edge  of  the 

boundary  layer.  The  second  momentum  equation  Indicates  that  the 
change  in  pressure  across  the  boundary  layer  is  small. 


EnetRy_ 


+ 5  Ci  5 

+  w |-(i parV*»‘  +/*  3 


(7) 


where  the  average  total  enthalpy  is  defined  as  It-  -  ^  K-  +"^’ 

i 

and  the  enthalpy  of  each  species,  ,  includes  the  chemical 

enthalpy  h? 
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It  is  n>w  useful  to  introduce  eddy  transport  coefficients 
for  turbulent  flow  . 


It  is  Lo  be  noted  that  fluctuations  in  fij  are  due  to  temperature 
only  so  that  the  definition  of  equation  (9)  is  reasonable.  Prandtl 
and  Lewis  numbers  for  turbulent  flow  are  then  defined  as? 


(U) 


f  Ce  (  O  +  iV) 

K  ■*-  *<-f 


(12) 
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Because  both  the  molecular  and  macroscopic  transport  properties 
arise  from  similar  mechanisms  for  many  ».ises  it  seems  reasonable  to 
assume  that: 


Po 


rs  La.  -  I 

i 


( 1  3> 


Under  these  conditiors  it  can  be  assumed  that  enthalpy  and  the  atomic 
mass  frictions  are  linear  functions  of  the  velocity  ratio.  With 
these  assumptions,  both  the  continuity  equation  for  atomic  mass 
fractions  and  the  enerpy  equation  reduce  to: 


_  LyL.«-  - 


's. 


1 


(14) 


This  equation  is  identical  wilh  the  momentum  equation  except  for  rhe 
term: 


A 


ptA.  e*  ^ 

ft  ' 


Although  neglect  of  this  term  in  the  moraeutura  equation  will  probably 
Have  little  effect  on  the  velocity  distribution  in  turbulent  flow, 
the  assumption  of  similarity  between  total  enthalpy,  atomic  mass 
fractions  and  velocity  ratio  is  somewhat  questionable  in  regions  of 
high  pressure  gradient.  By  the  reasoning  of  Lees,^  the  term  is 
probably  negligible  if  the  body  surface  is  strongly  cooled  because 
the  increase  in  density  near  the  wall  compensates  for  the  decrease 
in  velocity;  therefore,  the  coefficient  ot  the  pressure  gradient  is 
small  threugnout  the  boundary  layer.  This  argument  should  be  viewed 
with  some  caution,  however,  because  other  combinations  of  terms 
could  be  devised  such  that  the  coefficient  of  the  pressure  gradient 
would  be  small  throughout  the  boundary  layer.  Under  these  conditions 
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tfct  tertis  retained  in  the  momentum  equation  would  rot  lead  to 
similarity  between  the  enthalpy,  mass  fractions  and  velocity  ratio. 
The  advantage  of  similarity  is  that  once  the  velocity  distribution 
is  obtained,  the  use  of  the  von  Karma  n  momentum  integral  equation 
alone  is  sufficient  to  estimate  the  skin  friction,  heat  transfer 
and  mass  transfer  at  the  surface.  Otherwise  the  coiresponuiug 
integral  equations  fur  “neigy  ,-nd  species  continuity  :-iist  also  bo 
solved . 
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j.l  INTERACTION  OF  DIFFUSION  AND  CHEMICAL  REACTIONS 

In  a  chemical  system  in  which  diffusion  plays  an  important 
role,  use  of  the  continuity  equations  for  atomic  species  is  conven¬ 
ient  for  estimating  mass  transfei  rates  H-wever,  because  these 
equations  are  formed  free  linear  sums  of  the  molecular  species 
equations  they  usually  are  not  sufficient  to  completely  replace  the 
molecular  species  continuity  equations.  In  other  words,  there  are 
usually  less  atomic  than  molecular  species.  A  discussion  of  the 
number  of  additional  relationships  required  is  given  In  Appendix  B. 
These  additional  relationships  are  obtained  by  making  some  assump¬ 
tions  about  the  rates  of  production,  As  pointed  out  in 

reference  l,  the  specific  reaction  rate  coefficients  are  usually 
not  known  to  better  than  an  order  of  magnitude.  Therefore,  the 
study  of  limiting  cases  through  utilization  of  simplified  models  of 
the  actual  chemical  kinetics  is  more  profitable. 

In  the  first  place,  it  is  assumed  that  some  possible 
reactions  do  not  occur  at  all.  This  entirely  eliminates  some  of  the 
possible  i»o!ccular  species  from  further  consideration.  In  addition, 
two  opposite  limiting  assumptions  can  be  made  about  the  rates  of 
production,  Ur^  .  It  might  be  assumed  tb«t  the  specific  reaction 
rate  coefFlcients  are  very  small  so  that  reactions  car  occur  only  at 
the  solid  boundary  or  c  -tride  the  boundary  layer.  Under  these  condi¬ 
tions  the  molecular  mass  fractions  can  be  obtained  as  linear  functions 
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of  the  velocity  ratio  because  the  troublesome  U?i  terms  can  be 
neglected  In  equation  (1).  In  most  situations  involving  combustion, 
however,  the  high  temperatures  in  the  boundary  layer  probably  lead 
to  High  tale  coefficients.  The  opposite  limiting  case  is  to  assume 
that  these  rate  coefficients  are  very  large  compared  to  the  produc¬ 
tion  rates.  Under  these  conditions  the  perturbation  from  equilibrium 
caused  by  diffusion  and  convection  ot  species  is  small.  Therefore, 
the  mixture  ui  gases  at  any  point  in  the  boundary  layer  can  be 
determined  from  a  knowledge  of  the  atomic  mass  fractions,  T***,  and 
through  use  of  expressions  for  the  equilibrium  constants  of  the 
various  reactions.  This  procedure  will  be  adopter!  for  the  remainder 
of  this  report.  If  the  disturbance  due  to  diffusion  is  of  some 
importance  it  still  may  be  convenient  to  use  this  approach  for  Lhe 
first  iteration  in  the  computation  procedure. 


3.2  CHEMICAL  MOfrfcL 

Although  a  more  general  treatment  Is  possible,  it  is  now 
assumed  for  simplicity  that  the  surface  material  1$  a  sublimator, 
i.e..  It  passes  into  the  gaseous  state  without  going  through  a 
liquid  phase.  In  addition,  it  is  assumed  that  it  reacts  wi»h  oxygen 
to  form  two  possible  products  of  combustion,  P,  and  Px  ,  but  it 
does  not  combine  with  nitrogen.  The  possibility  that  oxygen  and 
nitrogen  may  dissociate  or  recombine  is  also  Included,  but  formation 
of  other  species  such  as  HO  and  ionized  molecules  is  neglected  It 
was  shown  in  reference  1  that  these  latter  assumptions  are  adequate 
for  air  over  a  wide  range  nf  temperature .  A  set  of  possible 
reactions  in  the  gas  phase  are  then: 


V0  0  f-  v^.  4- 


V  ’  P  4-V”  P 
,p,  r4  -t-Vp  r. 


(15) 


0zt;  zo 


(16) 


N  ,  ^  2  IN 


(17) 
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where  represents  the  .urface  material  (fuel).  II  i.i  furthtr 
assumed  that  the  fuel  dues  not  contain  nitrogen,  hi  t  may  conti 
oxygen  in  audit  i  in  to  two  atomic  species  -f,  and  .  The 

equilibrium  constants  corresponding  to  the  above  reactions  arc 
given  by: 


Vj'-V;' 


Tr 


i--  ^,0,4 


(18) 


(19) 


(20) 


The  partial  pressures  Pt  are  related  to  the  mass  fractions  K'L 
through : 

P  -  £i_  MP 

-  Mi 

where  M  is  the  molecular  weight  of  the  mixture  and  P  is  the 
pressure  for  the  mixture.  The  equilibrium  constants  lfj»„  are 
usually  i-nown  rather  accurately  as  functions  of  temperature. 


) . 3  MASS  TRANSFER  AT  SURFACE 

The  boundary  condition  at  the  surface  for  mars  transfer 
can  be  stated  in  terms  of  conservation  of  atomic  species: 


-  11  - 
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✓ 

•  c 

vhcre  sy*  *j  Is  the  net  flux  <'r  surface  material,  while  \jj  is 

the  net  flux  of  molecular  species  i  at  the  surface.  For  each 
reaction,  conservation  of  atomic  species  leads  to: 

A  =.  O  (23) 

U 

where  is  the  number  <•!  atoms  of  atomic  species  -w  per  molecule 

of  molecular  species  t  .  It  is  tctuLcd  to  through 

s  M***/>ML  ■  Simultaneous  solution  of  the  sets  of 

equations  given  by  equations  (22)  and  (23)  has  among  the  results  the 
following  expression  of  surface  mass  flow  for  the  particular  chemical 
system  under  consideration: 


This  equation  states  that  the  net  flux  of  surface  material  is  equal 
to  the  sum  of  net  fluxes  due  to  the  material  which  sublimes  and  that 
which  enters  into  chemical  reaction  with  oxygen  at  the  surface. 

The  net  flux  rates  consist  ot  the  difference  between  the 
flux  of  the  species  issuing  from  the  surface  and  that  which  reacts 
with  the  surface.  In  other  words: 

Kw  —  A*- 1  f  ~^L  ^ib  (25) 

*  .  •  , 
where  is  the  flux  of  species  L  from  the  surface,  is 

the  flux  striking  the  surface,  and  <TW  is  the  fraction  of  this  hack 
flux  which  reacts.  In  the  case  of  the  fuel  <J>  is  commonly  known 
as  the  sticking  factor  or  the  accommodation  coefficient.  The  back 
flux  can  be  calculated  easily  from  elementary  kinetic  theory^). 
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ub  — 


2  Tl*  G  ~T"mj 

m; 


T!k-  outward  rates  can  .i!s<  be  estimated,  but  they  arc 
better  treated  separately.  In  the  case  of  the  fuel,  the  outward 
rate  is  due  In  the  sub  1  in.a t  ion  process.  Vne  mass  fljw  rate  should 
depend  principally  upon  the  solid  surface  temperature  and  reduce  to 
no  net  mass  flow  at  equilibrium.  Therefore: 


’++  " 


-  %  K? 


1)  ?tM?o  T.J 


Here,  the  quantity  is  the  equilibrium  constant  for  sublima* 

tion.  An  expression  of  similar  form  could  be  deduced  from  the 
physics  of  the  solid  surface.  The  outward  rate  for  oxygen  has  a 
somewhat  more  complex  cause.  In  this  case,  some  of  the  products 
produced  by  reaction  between  oxygen  and  the  surface  may  encounter 
conditions  in  a  microscopic  layer  near  the  surface  such  that  they 
dissociate  again  and  precipitate  fuel  species  on  the  surface.  The 
result  is  that  oxygen  is  released  into  the  boundary  layer.  Clearly, 
such  a  process  must  depend  not  only  on  the  wall  temperature,  but 
also  on  the  presence  of  the  products.  In  order  to  reduce  to  no  net 
flow  at  equilibrium,  the  expression  for  atomic  oxveen  is: 


yfai 

Y  ko>,^ 

Nr*  v*# 


Similar  reasoning  for  molecular  oxygen  leads  to: 

r*  ^  n  z. 


-  'Si  1 


~q  ft  Po~r^7 

Mo, 
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Suhst  i  tin  I  op  ol'  equation*  ('.*>)  through  (2*J)  ii  equation  (26)  leads 
to  an  expression  [or  the  net  mass  flow  1 rom  the  surface  in  terms  of 
wall  temperature  and  the  partial  pressures  or  mass  fractions  of  the 
various  species.  Since  eqtii  I  ihrium  is  assumed  for  the  gas  phase 
reactions,  the  mass  fractions  are  related  through  equations.  (18), 
(lii),  and  (20).  Additional  relationships  are  obtained  frr.n  solutions 
of  the  boundaiv  layer  aqua tions ■ 


3.6  ATOMIC  MASS  FRACTIONS  AT  SUKFACE 

If  the  Lewis  number  is  unity  and  the  atomic  mass  fractions 
are  linear  functions  ot  the  velocity  ratio,  then  equation  (22)  leads 
te  the  following  set  of  expressions  for  atomic  mass  fractions  at  the 
surface : 


lTft 


(30) 


where 


?au-*  c+e/^ 


(30 


The  quantity  B  ,  which  represents  the  ratio  of  mass  injection  to 
skin  friction  or  heat  transfer  coefficients,  often  appears  in  combus¬ 
tion  literature.  A  further  relationship  between  skin  friction  and 
&  is  obtained  from  a  solution  of  the  momentum  equation  which  will 
be  discussed  in  a  later  secLion  for  turbulent  flow.  The  system  of 
equations  represented  by  equation  (30)  can  be  put  in  a  more  conven¬ 
ient  form  for  the  particular  chemicel  system  under  consideration. 
Simultaneous  solution  of  equations  (30)  and  (23)  leads  to  the 
following  equations: 


K”*j  n/ 


TV*  a. 

i  -eS 


(32) 
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K 


».  u; 


where 


V 

l  oe. 

Vo  Mu  1 

_“»S 

..  S-B* 

+-  ’M 

i  t-a 

V,' 

,  V,'M, 

v 'Mo 

1  ce 

V 


i-8 


( 33) 


(3A) 


The  above  equations  represent.  Sour  relations  among  the 
eight  unknowns:  KV,W^  KWlWj  ^  w/rffWj  Kew  j  *0**,  .and  6 

Three  additional  relations  are  supplied  by  the  equilibrium  relations 
of  equations  (18),  (19),  and  (20).  If  pe‘44-ci c/i.  ‘s  known,  one 
additional  relation  is  supplied  by  equations  (2A)  through  (29). 
Therefore  ail  eight  unknowns  can  probably  be  found  as  functions  of 
pressure,  wall  temperature,  and  the  quantity  ft*-t  C+i/a.  ■  Such 
relations  apply  to  either  laminar  or  turbulent  flow.  As  mentioned 
previously,  the  connection  between  and  the  other  variables  for 

turbulent  flow  will  be  discusced  in  a  iater  section. 


3.5  HEAT  TRANSFER 

The  heat  absorbed  by  or  radiated  away  from  the  surface, 

,  is  obtained  from  a  heat  balance: 

(35) 

The  heat  flux  is  due  to  conduction  out  of  the  boundary  layer  less 
the  difference  between  heat  convected  by  all  individual  species  and 
that  convected  from  the  solid  whose  surface  enthalpy  is  h+s  .  If 
Se  i  (Vc  -=  l  and  the  total  enthalpy,  ,  is  a  linear  function  of  the 

velocity  ratio,  then  it  can  be  shown  that: 
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"  kw '  B  (hs^,  ^^3  (36) 

P.y  making  use  of  expressions  (323.  (33),  and  (34),  4?^  can  he  put 
in  the  following  more  illustrative  form:  9 


|-6.-  ft^ciaA  jjn«.  +  ^e.Qp0w~  ^<w}(37) 


“'Ri’b-o  '  ^  1  y  ■  J|"o  W  *■  fow  +~  (■-£**) 

%  Pa. 

^ew  -  T*«t  (^M-v  hNw  4- 


+-T«t  (€„w  Uow  ♦(l-towH0J 

A*w  -  —  tn+s 

<*w  -  jj  ~  ^  fow+fitwVTse  ] 

60W  =  ^Ow/  (K®w+-Ko»w')  £-»«,-  ^*Jw/  (tiva’^wj 

The  quantity  Qpow  is  the  heat  ot  combustion  at  wall  conditions, 
h  is  the  enthalpy  of  air  at  wall  conditions  of  temperature  and 
dissociation,  and  Xf-vy  is  the  beat  of  sublimation.  When  there  is 
no  combustion  or  sublimation,  both  and  o^io  are  zero,  So  that 
equations  (36)  and  (37)  reduce  to  the  ordinary  results  for  compres¬ 
sible  flow  with  Prandtl  number  unity.  If  combustion  takes  place 
t^vu  approaches  unity  because  most  of  the  oxygen  is  used  up - 
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A.  I  VELOCITY  DISTRIBUTION 

Inspection  of  the  momentum  equation  given  by  equation  (5) 
shows  that  In  the  sublayer  the  shear  must  increase  from  the  surface 
value  due  to  the  etfect  of  mass  injection  ot  the  surface.  The 
Increase  in  shear  near  the  wall  can  be  estimated  byt®); 

'X  -  T*,  (  I  T  B 2-T  (38) 


whe  re 

-v  _  rr  k*S  _ 
L  ~  ^ 

(p.AT)  t*-' 

a  -  Km.  M  Me. 

« 

2  /M.  vij 

D  -w 

-  a  / 

2  *  /v-s. 
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The  quantity  \L  ,  which  has  appeared  in  previous  relationships,  is 
a  measure  of  the  ratio  of  mass  injection  to  skin  friction  or  heat 
transfer.  The  exponent  is  introduced  in  order  to  trace  the 
influence  of  the  shear  distribution  on  skin  friction.  Its  value  is 
unity.  This  shear  distribution  should  be  viewed  with  the  same 
qua  1  i f ications  as  apply  ter  the  constant  shear  assumption  usually 
made  for  turbulent  boundary  layers  with  no  blowing  The  shear 
should  approach  zero  at  the  outer  edge  o!  the  boundary  layer,  but 
the  influence  of  lire  vehreity  distribution  in  this  region  on  the 
skin  friction  may  he  srna 1 1 . 

In  the  fully  turbulent  region  near  the  wall,  it  is  assumed 
that,  following  Van  Priest^); 

✓ 

't'  ~  ~  f  AT*  *  (39) 


Introduction  of  the  Prandtl  mixing  length,  Q  ,  leads  to; 


/V  -  /.  2  j  bjX  1 

f  *  '  =>^1  ^ 


(sO) 


Either  the  assumption  of  Prajidtl ,  i  :K*,,  or  the  von  Kjrman  similar¬ 
ity  law,  JL  -s  K  may  now  be  introduced.  It 

is  assumed  that  these  formulations,  originally  developed  for  incom¬ 
pressible  flow  with  little  yL  variation  of  fluid  properties,  apply 
to  the  present  case.  By  integration  of  equation  (40)  with  the 
Prandtl  assumption,  the  velocity  distribution  is  found  to  he  given 
implicitly  by; 


V 


K  \Aq_ 


e_ 


(*i> 
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In  both  cases  the  arbitrary  functions  of  X  which  result  from  inte¬ 
gration  have  been  chosen  so  that  the  velocity  distribution  reduces 
to  that  of  incompressible  flow  with  no  Mowing  when  the  gas  proper¬ 
ties  and  shear  do  not  vary.  It  has  been  pointed  out  by  several 
writers'  ) 1  that  there  are  many  formulations  which  could  be  made 
to  reduce  to  the  incompressible  relations  depending  upon  where 
piuperties  in  the  tunctions  of  integration  are  evaluated.  The  above 
formulation  follows  that  of  Van  Driest'*^  who  based  these  proper¬ 
ties  on  wall  conditions  be-ause  of  rhe  strong  influence  of  the  wall 
in  the  sublayer.  .Justification  can  only  come  from  a  comparison  with 
experiments.  Van  Priest’s  results  for  compressible  flow  on  flat 
plates  have  given  rather  good  agreement  with  experiments  . 

The  change  in  variable  from  AJ  to  g  .at  constant  x.  , 
for  use  in  the  momentum  integral  is  given  by: 
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Vw- 

tv'  ^e. 


(ay-Kr)  r  J5 _ .  “j  ' 


(43) 


rfhero  is  zero  it  tho  von  iCarnkin  similarity  law  is  used,  and 

unity  if  Jt  —  K*^. 


4  2  SOLUTION  OF  MOMENTUM  INTEGRAL  FOR  SKIN  FRICTION 

The  von  Karman  momentum  integral,  modified  to  include  the 
effect  of  mass,  injection,  can  be  obtained  either  by  integration  of 
the  momentum  equation  across  the  boundary  lavei ,  neglecting  double 
and  triple  correlations,  nr  by  consideration  of  momentum  conserva¬ 
tion  for  an  element  ot  lluid.  The  result  is: 

<*<L  ^  /  f  *  /?°*(l  ~ 

+  t  UJ  (i  Tti)  (44) 


where  as  betor?  K.  is  unity  for  bodies  of  revolution  and  zero  for 
two-dimensional  bodies.  It  is  consistent  with  previous  assumptions 
to  neglect  the  first  term  on  the  right  hand  sfde  of  equation  (44). 

In  order  to  change  the  variable  of  integration  to  the  velocity  ratio, 
Z  ,  equation  (43)  is  substituted  In  equation  (44).  The  result  ran 
be  put  in  the  form: 


P  -  az  —  /a\  \\  - 

rr*  ct*  ^  j. 


(45) 
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where 


-K(- 

D  -  <? 

(  $ /  r  9  _ 1  (%. 

J  '  J  /?-  t-Tu^yM  *(-*)«.  d* 


2. 

The -quant  l tv  ft-  which  up  .nears  in  equation  (66)  is  inversely  pro¬ 
per  t  !«•::« !  to  the  shin  frit. ion  coefficient  based  on  wall  conditions. 
Since  lX3"  l:;  an  unknown  t unction  of  X  ,  equation  (65)  cannot  be 
integrated  as  it  stands.  However,  at  a  sufficient  distance  from 
the  stagnation  point,  £*"  should  he  large.  This  permits  certain 
approx inut lions  to  be  made.  Change  of  variable  from  2  to  makes 
It  possible  to  express  the  integral  J  as: 


where 


/ 


(1  +  82)“] 


(66) 


Successive  integration  ol  w)  by  parts  produces: 


g- 

a. 


O.  &  09  r  _  /n  (<*.  \ 


0\ 


B.  d  s 


] 


*<L 


(6  7) 


For  large  ft  with  (4^,  of  order  unity  the  dominant  term  yields: 

ft^i  \  r  o. .  ,  ~^2_ 

J-  jrt  e  (its)  [-ftT 


(68) 
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In  addition,  1 1  tan  be  -i  luiw.i  Lh.il  11  r.  (Js'{  is  lari'..-: 


U 


'J 


•:/  -V 


$ 


a  % 

l  vj 


j 


>;q.ia:!.i::  tin  now  be  inlejra  ted  approximately  by  making  use  of 

equations  («S)  and  (a1!)  Tile  icsul!  ni vis  a  law  for  variation  of 
turbulent  sl.in  friction: 


where 


Re  -  o 

A/  I  - 


a 

/  -2- 


^e2/<w,  />o^ 


Jt  has  been  assumed  that  the  influence  of  wall  temperature  and  mass 
fraction  variations  in  the  X-  direction  is  negligible.  If,  in 
addition,  it  is  assumed  that  the  variation  of  life  is  not  signifi¬ 
cant  in  determining  the  skin  friction,  the  effective  Reynolds  number 
becomes : 
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This  Reynolds  number  reduces  to  the  ordinary  Reynolds  number  based 
on  wall  conditions  or  a  Mat  plate.  for  a  conical  geometry,  where 
At,  —  C  X  ,  it  becomes  half  the  ordinary  Reynolds  number  consistent 
with  a  derivation  by  Van  Driest^'')  ior  the  conical  geometry .  The 
quantity  <0  is  unity  for  incompressible  flow  with  no  mass  injection, 
for  compressible  flow  witii  constant  specifii  heat  and  no  mass  injec¬ 
tion,  <p  reduces  to  Van  Driest's  result  when  the  enthalpy  velocity 
variation  is  introduced  In  the  -  ase  of  combustion,  the  evaluation 
oi  is  much  mote  complicated  due  to  the  influence  of  variations 

in  the  gas  mixture  through  the  boundary  layer  which  affect  the 
relationships  between  density,  enthalpy  and  veloci  ty  ratio.  A 
method  for  obtaining  <p  is  discussed  in  Section  *.3  and  calculations 
oi  0  for  the  case  of  graphite  combustion  are  described  in 
Section  5.3. 

Equation  (50)  shows  clearlv  the  inMueute  ol  shear  varia¬ 
tion  and  mixing  length  absorptions  on  the  skin  friction  through  the 
exponents  ot  and  f3  As  Is  well  known,  both  mixing  length  formu¬ 
lations  lead  to  the  same  result  for  incompressible  flow  with  .to  mass 
injection  Equation  (50)  has  the  same  fortr.  ,ts  the  incompressible 
von  Karman-Schotnherr  law  for  skin  friction.  Two  arbitrary  constants, 
K”  and  0  ,  appear  which  must  He  obtained  from  experiments.  If 
these  constants  are  assumed  to  be  the  same  as  those  obtained  for  the 
incompressible  case,  then  equation  (50)  can  be  expressed  as: 


I.  ?0  + 


(us) 


■Ki 


(52) 


It  would  be  truly  fortunate  if  these  constants  do  not  have  to  be 
adjusted 
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It.  should  he  pointed  on*,  th.it  because  of  the  ass  ..mpt  i  on  of 
large  <s*.  Lite  skin  filclion  iaw  does  not  apply  at  the  stagnation 
point  of  lilunt  bodies.  In  this  ronton,  the  oht>ve  formulation  leads 
to  finite  shear  stress  and  infinite  heat  transfer,  both  of  which  arc 
phy ::  it.:  1  !  y  unsound.  However,  since  the  Mow  will  probably  not 
become  turbulent,  until  at  least  some  small  distance  from  the  st.ign/i 
lion  point,  the  behavior  in  this  region  is  not  as  important  as  in 
the  case  of  laminar  flow. 


'>.!  r.tTKCT  OK  CHEMICAL  STRUCTURE  OP  BOUNDARY  LAYER  ON  INTEGRAL  (f) 

All  of  the  results  obtained  so  far,  except  for  evaluation 
of  the  integral  0  ,  have  required  only  evaluation  of  fluid  proper¬ 
ties  at  the  body  surface  or  at  the  edge  of  the  boundary  layer.  In 
the  case  of  laminar  flow  with  assumptions  similar  to  those  made  in 
the  present  investigation,  no  detailed  knowledge  of  the  chemical 
structure  ot  the  boundary  layer  i3  required  in  order  to  estimate 
beat  transfer,  mass  transfer  and  skin  friction.  However,  with  tur¬ 
bulent  flow,  it  is  not  possible  to  make  the  same  transformations  of 
the  boundary  layer  equations,  because  the  transport  coefficients 
are  not  properties  of  the  fluid  but  depend  on  the  flow  field.  For¬ 
tunately,  with  the  assumptions  which  have  been  made,  the  only 
difficult  quantity  to  compute  is  the  integral  0  ,  Its  evaluation 
requires  a  knowledge  ol  the  density  variation  through  the  boundary 
1 ayor . 


When  diffusion  and  chemical  reactions  occur  in  the 
boundary  layer,  the  density  varies  not  only  because  of  temperature 
variations /^but  also  because  of  changes  in  chemical  composition. 
The  temperature  itself  is  related  to  the  enthalpy,  which  is  known 
in  terms  of  the  velocity  ratio  In  a  rather  complicated  implicit 
tashion.  Not  only  is  Q)  more  complicated  Co  compute,  but  in 
addition,  if  depends  upon  more  independent  parameters. 

The  last  three  parameters  must  be  added  in  order  to  determine  the 
composition  throughout  the  boundary  layer.  Apparently,  a  simple 
graphical  representation  of  0  is  not  possible.  However,  the 
Integrand  ot  0  can  be  determined  if  the  enthalpy,  pressure,  and 
gas  composition  are  known. 
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Fortunately,  with  the  assumpt  ior:s  of  fast  reaction  rotes 
In  the  >ms  phase  and  LC  -  I  ,  the  gas  composition  can  be  determined 
in  terms  of  n  single  parameter,  £  .  This  can  be  seen  through  use  of 
equation  (JO)  for  the  atomic  mass  fractions  at  the  surface  and  use 
of  a  linear  relationship  with  the  velocity  ratio,  introduction  of 
the  parameter  5  <  defined  below,  then  yields  the  following  exptes* 
slon  for  a  tor**’  -:i  ,s  fractions  at  any  point: 


where 


f  -  ('-2)4s 


(53) 


In  order  to  actually  compute  the  density  variation,  the  molecular 
mass  fractions  must  be  found.  These  are  obtained  by  simultaneous 
solution  of  equations  (53)  and  (23).  For  the  chemical  sysLem  under 
consideration  in  this  report  the  results  are: 


VTj.  =  J  —  +  %"t-  (If,  r*K)  (54) 

Vo  M  o 


if;.  *  Hi  f 


V0'Mo 
r'w 


j Jot 


(55) 


(56) 


In  addition  to  these  expressions,  the  equilibrium  constants  given  by 
equations  (18),  (19),  and  (20)  serve  to  determine  the  system. 
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T'ic  relationship  between  f  and  the  velocity  ratio  in  the 
bounda ry  layer  depend:.’  on  the  wall  conditions  of  mass  ttansfer  whicn 
determine  S  .  When  tut-  molecular  composition  is  obtained  at  each 
location  in  the  boundary  layer,  then  the  integrand  of  0  can  be 
evaluated  in  terms  of  enthalpy  and  pressure  and  0  computed 
numerically  from: 


I 


(57) 


It  should  he  noted  that  the  enthalpy  is  related  to  the  wall  and  free 
stream  boundary  conditions  and  the  velocity  ratio,  2  ,  through: 

-  (h<,c 2-  (58) 


As  stated  above: 


f 

L  P  h-ST 


i  (p  k  i) 


(59) 


The  numerical  computations  arc  Tided  by  determining  the  integrand 
*f“  as  a  function  of  pressure  and  enthalpy  at  various  fixed  values 
of  $  It  can  be  seen  from  the  definition  of  j  that  it  is  zero 
at  the  outer  edge  of  the  boundary  layer  and  increases  toward  the 
surface  when  there  is  mass  injection-  The  surface  value  of  the 
composition  parameter,  =8/('*8j  ,  can  be  obtained  in  terms  of 
the  mass  transfer  conditions  at  the  surface.  Under  assumptions 
previously  made,  its  value  depends  upon  wall  temperature,  pressure, 
und  the  parameter  ft  Thus  evaluation  of  0  requires 

simultaneous  solution  of  equations  (52)  and  (57). 
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Ii  addition  to  fyj  another  impo’-Uint  value  of  f  Ls  that 
which  occurs  at  the  place  where  the  reactants  ore  in  stoichiometric 
proportion.  Kquntion  (54)  shows  that  this  value  of  J  la  given  by: 


I  +  6* 


(60) 


This  value  of  ^  occurs  in  the  interior  of  the  combustion  zone,  a 
region  where  the  major  portion  of  the  combustion  reaction  takes 
place.  For  most  combustion  reactions  the  equilibrium  constant  given 
by  equation  (18)  is  very  small  for  the  environmental  conditions 
which  occur.  Therefore,  in  the  combustion  zone  the  mass  fractions 
of  the  reactants  are  small.  Under  these  conditions  the  model  of 
reference  l  is  well  approximated.  The  boundary  layer  is  divided 
into  two  regions  by  n  thin  combustion  zone.  On  one  side  of  this 
zone  oxygen  is  mixed  with  products  and  inerts  while  on  the  other 
aide  fuel  is  mixed  with  products  and  inerts.  An  estimate  of  the 
thickness  of  the  combustion  zone  for  the  case  where  combustion  takes 
place  in  the  gas  phase  is  given  in  Appendix  C.  If  wall  conditions 
arc  such  that  the  reactants  are  in  stoichiometric  proportion  at  the 
surface,  then  5u/ -  ¥*  •  At  lower  values  of  the  excess  of 

oxygen  increases  at  the  wall  until  finally  the  combustion  reaction 
docs  net  take  place  at  all  If  ^:o  . 


s  ten  on  i 

COMBUSTION  OF  GRAPHITE 


5.1  CHEMICAL  CONSIDERATIONS 

As  in  reference  i  the  combustion  ol  graphite  is  investi¬ 
gated  numerically  because  a  great  deal  of  data  is  available  on  this 
substance.  The  combustion  process  is  assumed  to  be  adequately 
described  by  the  single  reaction  step: 


C  +  O  X.  CO  (61) 

The  formation  of  CO2  which  will  only  occur  in  appreciable  quantities 
at  low  temperatures  is  neglected  In  conformity  with  previous 
assumptions  reaction  with  nitrogen  which  requires  extremely  high 
temperatures  is  also  neglected.  In  addition,  the  form  of  the  carbon 
molecule  in  the  gaseous  state  is  assumed  to  be  monatomic,  although 
recent  experiments,  (12)  (13)  indicate  the  possibility  of  C3  forma¬ 
tion  during  sublimation.  Since  this  question  is  not  entirely 
settled,  the  simpler  configuration  is  assumed  for  purposes  of 
illustration.  All  of  the  above  complications  could  be  included 
with  slight  modilication  of  the  analysis . 

The  basic  data  listed  in  Table  1  of  reference  1  are 
utilized  tor  the  present  calculations  Additional  data  are  required 
because  of  the  somewhat  more  general  treatment  of  the  combustion 
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process  in  ine  present  analysis.  Equilibrium  constants  fitted 
empirically  to  data  of  reference  (14)  are  as  follows: 


P  -  Kp,  ^  l,3iUx 


S  „ 

>o  e. 


_  /.Si  3  7Z  »./0 


(62) 


Pc. 


^  f 


7  370Q7x/»7S~ 

0  i  72  x  <0  £  T 


(63) 


The  equilibrium  constants  are  in  atmospheres  while  the  temperature 
is  in  degrees  Rankine. 

The  sticking  factor  or  accommodation  coefficient,  *3T  ,  for 
graphite  sublimation  is  taken  as  0.3.  This  fits  the  data  of 
references  (15)  ar.d  (16)  over  a  moderate  temperature  range  (up  to 
about  3000°K)  assuming  that  the  vapor  is  monatomic  carbon.  In 
reality  the  accommodation  coefficient  probably  varies  with  tempera¬ 
ture  but  the  single  constant  value  is  consistent  with  the  present 
state  of  knowledge.  The  value  of  <j  for  monatomic  oxygen  is 
assumed  to  he  unity  while  that  for  diatomic  oxygen  is  chosen  at  the 
two  extremes,  zero  end  unity,  for  purposes  of  illustration, 
specially  designed  experiments  are  required  in  order  to  determine 
these  coefficients. 


5.2  GENERAL  SURFACE  MASS  AND  HEAT  TRANSFER 

Calculations  have  been  made  for  graphite  by  means  of 
several  IBM  650  programs  using  the  analytical  results  obtained  in 
Section  3.  Since  nonequilibrium  surface  conditions  have  been 
considered,  it  is  necessary  to  know  a  flow  parameter  in  addition  to 
surface  temperature  and  pressure-  ir.  order  to  specify  the  burning 
rate.  This  parameter,  P«w«  ,  is  a  Erasure  of  the  transport 

of  properties  to  or  from  the  wall.  For  a  particular  body  geometry 
with  given  local  free  stream  conditions  as  well  as  wall  temperature 
and  pressure  this  parameter  car.  be  obtained  through  use  of  the  skin 
friction  law  of  equation  (52)  for  turbulent  flow  in  conjunction  with 


ASX  PUBLICATION'  NO.  11-166 


COMBUSTION'  OF  GRAPHITE 


✓ 

equation!,  (18)  through  (20),  (24)  througli  (2V)  ,  and  (32)  through 
(34)  ,Ior  the  mass  fractions  at  the  surface  and  1$.  In  this  section 
Lhc  quantity  is  treated  ns  an  independent  variable  so 

that  Hie  insults  apply  to  any  geometrical  configuration  in  laminar 
or  turbulent  flow.  Ail  of  the  above  named  equations  except 
equation  02)  are  utilized.  In  Section  5.5  calculations  are 
described  for  the  particular  case  of  a  hemisphere  in  turbulent  flew. 

In  Figure  1  the  mass  injection  parameter,  B,  is  presented 
as  a  function  of  wall  temperature,  flow  pn<*nr.u»ter,  and  accommodation 
coefficient  for  a  fixed  pressure  of  one  atmosphere.  The  entire 
region  below  B  *  B*  represents  conditions  under  which  the  combustion 
process  is  controlled  by  chemical  reuction  al  t!.e  surface.  In  this 
region  an  appreciable  quantity  of  oxygon  is  not  consumed  by  the 
reaction.  At  low  temperatures  the  value  of  the  acecumnodatlnn  coeffi¬ 
cient  for  molecular  oxygen  has  considerable  influence  on  the  burning 
rate  while  at.  higher  temperatures,  when  most  of  the  oxygen  becomes 
dissociated,  it  has  little  effect.  All  of  the  curves  swing  almost 
vertically  upward  at  still  higher  temperatures  due  to  the  sublima¬ 
tion  of  gaseous  carbon.  The  sublimation  rate  eventually  becomes 
sufficiently  large  to  drive  the  combustion  zone  away  from  the  surface 
so  that  B  B*  Under  these  conditions  the  diffusion  of  oxygen 
toward  the  wall  is  not  sufficient  to  consume  the  flux  of  carbon 
comLng  oft  due  to  sublimation. 

The  lower  values  of  the  flow  parameter  result  in  higher 
values  of  B.  When  the  flow  parameter  becomes  sufficiently  low  it 
has  no  influence  at  all  so  that  equilibrium  is  approximated.  At 
high  values  of  the  flow  parameter  the  quantity  B  becomes  very  small. 
This  dependence  of  B  on  the  flow  parameter  has  a  significant  effect 
on  the  heat  transfer  rate  given  by  equation  (37).  The  factor 
exhibits  a  similar  behavior  as  shown  in  Figure  2.  The  net  effect  of 
these  factors  is  presented  in  Figure  3,  where  a  portion  of  the 
enthalpy  potential,  namely,  ,  is  presented 

as  a  function  of  t'nc  wall  temperature  and  flow  parameter  at  one 
atmosphere  pressure.  The  flow  parameter  actually  has  the  opposite 
effect  on  the  surface  mass  loss  as  shown  in  Figure  4.  Higher  values 
of  the  flow  parameter  result  in  higher  surface  mass  losses. 

The  effect  of  pressure  on  B,  and  enthalpy  potential  is 
shown  in  Figures  5,  6  and  7.  Increased  pressure  has  an  effect 
similar  to  that  of  decreased  flow  parameter.  Equilibrium  is  more 
nearly  approached  because  of  backscattering  toward  the  surface. 
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5 .3  EVALUATION  OF  THE  INTEGRAL  (jj  FOK  GRAPHITE  IN  AIR 

Ry  making  use  of  equations  (54)  through  (56)  in  addition 
to  equations  (18)  through  (20)  the  molecular  mass  fractions  can 
be  obtained  as  functions  of  tei-:oerature  and  pressuie  Cor  constant 
values  of  J  .  The  entha.py  of  the  mixture  can  be  substituted  Cor 
temperature  as  a  variable.  Therefore,  the  integrand,  /P(i'l)*  J  , 

which  occurs  in  the  expression  for  ^  given  by  equation  (57)  can  be 
found  as  a  function  of  the  enthalpy  of  Lhe  mixture  and  tne  pressure. 
Since  the  molecular  species  are  forced  to  include  C,  CO,  0,  Oj,  N, 
and  N2  the  enthalpy  of  the  mixture  will  not  necessarily  be  zero  at 
zero  temper' ture  because  of  the  heats  of  formation  of  C  and  CO.  A 
shift  of  this  datum  so  that  all  enthalpies  are  zero  at  zero  tempera- 
turemakes  it  possible  to  obtain  a  universal  curve  for  the  integrand 
of  (b  over  a  considerable  range  of  pressure,  enthalpy  and  £  .  The 
heat'of  formation  for  the  mixture  can  he  approximated  by  expressions 
linear  in  ^  .  When  f *  ,  Kt  =*  O  at  low  temperatures  so  that 

from  equations  (54)  and  (55): 

s*r  <«> 

If  f2$\  (ffo  4»V)~0  so  that: 


&  L<  c  -  r..  k:lf-Dvtf‘-s*3«> 


The  relevant  numbers  for  graphite  have  been  substituted  for  the  heats 
of  formation  and  the  results  of  calculations  for  the  integrand  using 
the  full  set  of  equations  (54)  through  (56)  and  (18)  through  (20) 
with  the  mixture  heat  of  formation  subtracted  out  have  been  plotted 
in  Figure  8.  Values  of  J  equal  to  zero  which  represents  air,  J  * 
which  represents  the  combustion  zone,  and  2.  J  *  which  represents  a 
fuel  rich  condition  were  chosen. 

A  fit  to  this  distribution  was  obtained  by  adding  a  correc¬ 
tion  term  to  that  which  results  for  air  treated  as  a  perfect  gas. 

The  expression  for  the  fit  is  as  follows: 
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IP(',U  H  ‘ 


•  Ir  C.  +8  X  'O  rC  H)  ‘ 


(H  = 


where  H  is  the  enthalpy  with  a  shifted  datum  given  In  BTH/lb.  A 
slignliy  better  fit  may  have  been  obtained  with  another  form  fcr 
the  correction  term,  but  the  one  chosen  allows  easy  analytical 
integration  of  0 


The  result  of  integration  is  as  follows: 


if  B  -£  B* 


_  f\  (MO.  7«S->  ,</*  )  7-0.  (*»‘l£d2)  \ 

fie..  ‘•i  '  1 1 — ; -  ir-z —  1 


s-:>:7iS*  /o'4  \  (fc.yv  -A-  a/4  I J  «^> 


1 1  ri  t’  v 


4  .  -  j  '  i  r ?-9  «/ o'*  fb*+4«V.]  )  J-vYs.h  ‘  br.  \ 

.4f/  Vx  wu.y4 

*  0,‘lii  *'•}  f  £  Cj_g  )_  +  b^(C.vY/t  ■> 

4  c.  —4*7  s 

+( 1  .f.w.wtuv.,.,; }  1 fs„- (Vi!  . 

q  1/JC.v.  'V7~=.  I 

•»  c.).r.o-*f  MS’ 

L  ^ '•I  arc,  v  i1  J  rp,a\ 
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Although  the  above  expressions  for  (p  appear  rather 
formidable  they  actually  represent  a  great  simplification  in  that 
they  are  analytic.  As  pointed  out  in  Section  4.3  0  is  a  function 
of  five  independent  parameters. 

It  should  be  noted  that  when  B  =  0  equation  (t>7 )  gives  a 
relationship  for  turbulent  flow  of  air  over  nonreacting  surfaces 
including  the  effects  of  dissociation.  At  low  temperatures  under 
these  conditions  the  results  reduce  to  those  of  Van  Driest. 
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3 lNVISCL'i  ll.O  i-'II'LD  Fuk  BLUNT  BODIES 

For  purposes  ol  illustration  heat  and  mess  transfer  on  a 
hemisphere  will  be  described  in  Section  5. '3.  The  hemisphere  is  a 
simple  geometry  representative  of  blunt  body  shapes.  A  detailed 
investigation  os  the  inviscid  flow  field  cn  blunt  bodies  is  beyond 
the  stupe  of  this  report  but  it  deserves  attention  because  turbulent 
heat  and  mass  transfer  are  more  sensitive  to  the  flow  field  than  the 
laminar  rates  A  review  and  extension  ol  this  problem  was  recently 
presented  by  Van  Dvke.^’*)  He  points  out  that  for  the  hemisphere, 
modi  tied  Newtonian  flow  as  proposed  by  Lees^*®)  gives  good  agreement 
wiLh  his  more  elaborate  numerical  calculations,  at  least  up  to 
about  fifty  or  sixty  degrees.  Lees^^)  also  presents  some  experi¬ 
mental  correlation  with  Newtonian  flow  for  the  hemisphere.  In  view 
of  these  comments  it  appears  that  use  of  modified  Newtonian  flow  on 
a  hemisphere  gives  an  adequate  description  of  the  pressure  distri¬ 
bution  fur  use  in  illustrating  the  heat  and  mass  flux  distribution. 
The  pressure  distribution  is  therefore  given  in  terms  of  angular 
distance  from  the  stagnation  point  by: 

P/p  ^  p/*>  h  6~  Vpre)  CuS  9 

,rTi  '  rTt~  J  (69) 


p.,  i 

The  pressure  ratio,  ‘  y'  F,  ,  across  the  normal  shock  is  affected 
very  little  by  dissociation.  However,  there  i»  a  strong  effect  on 
the  density  ratio,  SV-yy,  ,  which  is  required  in  nrder  to  deter¬ 
mine  the  heat  transfer.  Plots  of  the  variation  of  these  quantities 
with  free  stream  velocity  and  altitude  is  given  by  Feldman  tn 
reference  19-  The  pressure  and  density  serve  to  determine  the  state 
of  the  gas  at  the  stagnation  point.  In  order  to  determine  condi¬ 
tions  at  some  other  point  on  the  body  the  assumption  that  the  air 
entering  the  boundary  layer  is  at  approximately  the  same  entropy 
level  can  be  used  in  addition  to  equation  (69)  for  tire  pressure 
distribution.  Other  thermodynamic  variables  can  then  be  determined 
from  Moilier  charts.  However,  when  pressure  variations  are  moderate 
an  approximate  method  is  simpler.  This  method  consists  of  intro¬ 
ducing  a  polytropic  pressure  density  variation  ac  constant  entropy, 
l.e.,P/j^**“  =  const,  where  the  exponent  is  defined  by: 
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A  plot*  of  vs  free  stream  velocity  is  given  in  Figure  9  for  the 

pressure  limits  shewn.  The  value  oi  Vfi.  Is  only  of  importance  in 
regions  far  uw.»*  from  the  stagnation  point  because  at  small  distances 
the  density  is  approximately  constant. 

With  these  assumptions  the  local  mass  flow  and  kinetic 
energy  which  occur  in  the  heat  transfer  expression  are  given  by: 


?€  e<t 


UJ  ^  ^  L  fe)  J 


Another  required  quantity  is  the  effective  Reynolds  number 
for  use  in  obtaining  the  skin  friction  coefficient.  For  the 
von  Kantian  similarity  law,  with  (3=0,  this  parameter  for  a 
hemisphere  is  given  by: 


*  The  author  is  Indebted  to  Dr.  S.  Lemper t  of  Aeronutronic  Systems, 

T  *•>  —  f*>*»  eiiArttvfnn  nranV 
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In  arriving  at  this  relationship  variation*-  in  wall  temperature  ana 
r-iolecular  weight  along  the  body  have  been  neglected.  The  quantity 
is  the  free  stream  Reynolds  number  based  on  the  radius  of  the 
iiviiiiuphere .  The  parameter  R  is  plotted  as  a  function  of  Mach  number 
and  position  in  Figure  10-  If  the  wall  temperature  and  the  free 
stream  conditions  upstream  of  the  shock  are  known  the  effective 
Reynolds  mrahtr,  *  can  be  rounc1  at  any  boJy  location. 


3  -  5  HEAT  AND  MASS  TRANSFER  DISTRIBUTIONS  ON  \  HEMISPHERE 

Without  adequate  experimental  data  it  is  too  much  to 
expect  that  the  present  semi -empirical  analysis  will  give  accurate 
quantitative  predictions  of  heat  and  mass  transfer  However,  if 
good  judgment  has  been  used  in  choosing  the  important  chemical 
reactions  anti  making  physic^!  assumptions,  the  trend  of  these 
phenomena  may  be  predicted  by  making  use  uf  the  analysis. 

Forced  convection  heat  transfer  on  blunt  bodies  with 
nonreacring  surfaces  exposed  to  hypersonic  inviscld  flow  conditions 
has  received  considerable  attention  in  recent  years.  It  is  there¬ 
fore  interesting  to  determine  what  the  heat  and  mass  transfer 
distributions  may  be  under  similar  external  flow  conditions  when 
the  surface  reacts  with  the  environment-  The  combustion  of  a 
graphite  hemisphere  exposed  to  free  stream  conditions  corresponding 
to  Mach  10  flight  at  altitudes  of  IUU.00C  or  50,000  ft  has  been 
chosen  as  an  example.  In  addition  to  Mach  number  and  altitude,  the. 
free  stream  Reynolds  number  and  the  wall  temperature  must  be  given 
in  order  to  specify  the  beat  and  mass  transfer.  Two  Reynolds 
numbers  an  order  of  magnitude  different  have  been  chosen  tor  each 
altitude.  The  wall  temperature  was  varied  from  2000  to  6000°R . 

Some  assumptions  about  the  acconmodation  coefficients  must  also  be 
made.  These  were  chosen  as  <TC  «  0.3,  *  l,  and  =  0  for 

purposes  of  Illustration. 

In  order  to  carry  out  the  computations  the  inviscld  flow 
field  considerations  of  Section  5.4  were  combined  with  equations  (18) 
through  (20),  (24)  through  (29),  and  (32)  through  (34)  for  the  mass 
fractions  at  the  surface  and  B;  equation  (52)  for  the  turbulent  skin 
flicrlon  law;  cquatior.3  (67)  and  (68)  for  ;  and  equation  (37)  for 
the  heat  transfer.  Use  of  the  turbulent  skin  friction  law  eliminates 
the  flow  parameter,  o *  /  X-  ,  as  an  Independent  parameter 

for  the  particular  case  being  studied.  The  computations  were 
accomplished  by  means  of  an  IBM  650  program. 
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Figures  11  through  14  contain  distributions  of  heat  flux, 
mass  tiux,  B  arffi  at ,  respectively,  for  the  100,000  ft  altitude 
while  Figures  15  through  18  contain  corresponding  distributions  for 
the  50,000  ft  altitude.  Comparison  of  these  curves  shows  that  the 
level  of  both  heat  transfer  and  mass  flux  is  a  1  must  on  order  of 
magnitude  higher  at  the  lower  altitude  at  approximately  the  same 
Reynolds  number  and  surface  temperature.  This  is  due  to  the  strong 
influence  of  the  free  stream  density  in  turbulent  flow  which  is 
reflected  in  the  local  value  of  .  Because  the  skin  friction 

is  only  slightly  affected  by  variations  in  Reynolds  number,  the  heat 
and  mass  transfer  rates  vary  little  with  Reynolds  number  at  low 
temperature.  At  high  temperature  the  general  heat  and  mass  transfer 
plots  of  Figures  3  and  4  indicate  that  there  is  much  more  sensitivity 
to  the  flow  parameter  for  <r0x.  -  o.  Therefore,  variations  in 
Reynolds  number  have  a  moderate  influence  at  high  temperature  as 
shown  in  Figures  11  through  18. 

The  heat  flux  variations  between  6000  and  8000°R  shown  in 
Figure  11  for  the  100,000  ft  altitude  show  a  rather  confusing 
pattern  which  can  be  explained,  however,  by  a  comparison  with  the 
corresponding  conditions  Cor  B  and  o <  presented  in  Figures  13 

and  14.  The  confusion  in  this  temperature  range  is  caused  by  the 
rapid  occurrence  of  gas  phase  combustion  as  illustrated  by  Figure  13 
and  the  rapid  consumption  of  oxygen  at  the  surface  as  shown  by 
Figure  14  when  the  temperature  is  increased.  The  heat  flux  varia¬ 
tions  at  50,000  ft  display  a  more  regular  appearance  as  shown  in 
Figure  15  because  gas  phase  combustion  is  not  approached  in  the 
temperature  range  shown.  The  prevention  of  gas  phase  combustion  at 
che  1 ower  al titude  is  primarily  due  "o  the  effect  of  increased  flow 
parameter. 


At  low  temperatures  the  heat  flux  distributions  in  both 
Figures  11  and  1.5  exhibit  a  peak  value  in  the  vicinity  of  60  degrees 
on  the  hemisphere.  This  behavior  is  characteristic  of  turbulent 
heat  flux  distributions  or.  r.on-rcacting  blunt  bodies  at  hypersonic 
speeds.  Under  these  conditions  the  enthalpy  potential  is  only 
slightly  affected  by  surface  conditions  and  the  skin  friction  cueffi 
cient  varies  moderately  with  geometrical  location.  Thus,  the  flow 
parameter,  fa.  U  C+4./JL  ,  and  hence  the  heat  transfer  reflect  the 
variation  in  around  the  body. 

The  mass  flux  variations  shown  in  Figures  13  and  16  also 
exhibit  a  peak  flux  near  40  degrees  in  many  cases  in  spite  of  a 
different  shape  in  the  B  distributions  of  Figures  14  and  17. 
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Iti  Figure  !9  the  variation  of  the  integral  (fi  around 
the  hody  is  shown.  The  value  ,.f  (h  is  remarkably  close  to  unity, 
i he  •  iKompressitle  value-  This  is  Decause  at  high  temperatures  the 
integrand  of  <p  is  mainly  in  the  flat  portion  of  the  curve  shown 
in  Figure  8  throughout  the  boundary  layer. 
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APPENDIX  A 


lu'RBUIXNT  CONSERVATION  STATIONS 
FOR  REACT  INC.  fAS  MIXTURES 


The  turbulent  conservation  equations  have  been  stuuiuv  tor 
compressible  flow  In  references  (2)  and  (20).  A  similar  investiga¬ 
tion  is  carried  out  here  for  reacting  gas  mixtures  in  a  slightly 
different  manner. 

The  instantaneous  value  of  each  of  the  unknowns  is 
separated  into  a  mean  and  a  fluctuating  part  such  that: 


2 


-  ^  f-2  ‘ 


(A-l) 


where  may  be  any  one  or  any  couiblnalion  of  the  physical  variables. 
Introduction  of  expressions  such  as  equation  (A-l)  into  the  full 
instantaneous  conservation  equations  ior  reacting  gas  mixtures  and 
averaging  according  to  the  Reynolds  rules  tor  ti  uvetages  yields 
the  following  conservation  equations  which  inclu  correlations 
between  fluctuating  components: 


Continuity  for  Molecular  Species 

'B  *v-[(ff)K  MP?)V  *? (Bti, 


(A -2) 
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.■here  Tv  if  lilt-  muss  (riittiui,  the  diffusion  velocity  aiu 

Hr  is  the  mass  r-'te  of  productioi/  per  unit  volume  of  Lne  lt(i 
spot  ta.  Lt  is  assumed  that  diffusion  due  to  temperature  and 
pressure  gradients  is  negligible  compered  to  that  due  to  concentra¬ 
tion  gradients  and  that  all  dilfusior  coefficients  are  equal.  ! . •  r 
these  conditions  the  average  of  molocular  diffusion  terms  can  be 
expanded  to  g i vo : 


0  7^  -  p1  O'VK-  •-  (  pl>‘eP‘D  r  P  U')7lf‘CA*3) 


where  0  is  the  diffusion  cool  1"  icien  t .  Two  additional  continuity 
equations  cun  be  formed  from  lineat  sums  of  equation  (A-2). 

Overall  Continuity  Equation 

rV-(pf)  - 


(A-4) 


This  equation  is  obtained  by  suraming  equation  (A-2)  over  l  . 
Continuity  1  or  Atomic  Species 


l  -» '  fT  /  r*\  bj— i  * 


-o  (A-5) 


where 


t 


■  is  the  atomic  mess  fraction  of  species  ■4*t  Irrespective  of 
molecular  configure l ion  and  ■'l.tn.i  is  fbo  mass  of  atomic  species 
a-t  per  unit  mass  of  molecular  species  t  .  This  equation  is 
obtained  by  multiplying  each  expression  (A-2)  it*  and  summing 

over  t  .  The  assumption  uT  equal  diffusion  coefficients  leads  to 
the  tallowing  expansion  for  the  diffusion  terms  in  equation  (A-5): 


ASI  PUBLICATION  NO.  U-166 


VURRULEUT  CONSERVATION  E  NATIONS 


where  Is  the  coefficient  of  viscosity,  X  is  the  unit  tensor, 

q  is  the  rate  of  strain  tensor  and  p  is  the  pressure. 
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+  pf(S^‘ 

+  '?p'(fKV +£)]--  o 


where  In^  Is  the  total  enthalpy  including  chemical  enthalpy  and 
kinetic  energy,  h;.  is  the  enthalpy  of  ihe  dh  molecular  species 
including  chemical  energy,  and  Q>  is  the  flux  of  heat  by  conduction 
and  diffusion.  The  average  of  molecular  dissipation,  work  and  heat 
transport  is  given  by: 
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Equation  of  State 


where 


M 


l 


(A-1L) 


In  each  conservation  relation  given  by  equations  (A-2)  , 
(A-5)  .  (A-7)  and  (A-9)  the  first  term  inside  the  brackets  represents 
the  influence  of  convection,  while  the  second  term  represents 
transport  of  properties  by  turbulent  motion.  In  the  case  of  the 
momentum  equation  this  latter  term  represents  the  Reynolds  stress 
tensor . 


The  third  terms  whi<~h  involve  averages  of  molecular  trans¬ 
port  quantities  have  been  expanded  in  equations  (A-3) ,  (A-6) ,  (A-8) 
and  (A-10).  In  each  of  these  equations  the  first  terms  or  group  of 
terms  do  not  involve  any  fluctuating  components.  Near  a  body  surface 
in  the  laminar  sublayer  where  all  fluctuations  are  damped  out  due  tn 
the  influence  of  the  surface  and  property  gradients  are  large,  these 
first  terms  are  dominant  ones.  In  the  fully  turhulent  region,  how¬ 
ever,  the  turbulent  transport  terms  overshadow  all  components  of 
molecular  transport.  Presumably,  in  this  region  the  various  compon¬ 
ents  of  the  average  molecular  transport  terms  which  involve  correla¬ 
tions  between  fluctuating  transport  properties  or  fluctuating  gradi¬ 
ents  and  other  fluctuating  quantities  will  be  no  larger  in  magnitude 
than  the  largest  component  of  the  corresponding  term  involving  only 
average  transport  properties  and  gradients  in  equations  (A-3),  (A-G) , 
(A-8) ,  and  (A-10).  Since  it  is  to  be  expected  that  all  components 
of  the  turbulent  transport  terms  are  about  the  same  magnitude  each 
component  involving  correlations  between  molecular  transport 
fluctuations  will  be  overshadowee  by  the  turbulent  transport  terms 
In  the  fully  turbulent  region.  Presumably,  corresponding  deriva¬ 
tives  behave  the  same  way.  Therefore,  it  should  be  necessary  to 
retain  only  the  Tirst  term  or  group  of  terms  in  Ibe  expansions 
given  by  equations  (A-3),  (A-b) ,  (A-8)  and  (A-10),  except  possibly 
in  the  buffer  layer  between  the  sublayer  and  the  fully  turbulent 
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region.  No  aLlempl  will  be  made  to  account  tor  this  latter  compli- 
t.i  i  I  on  in  the  buffer  layer.  The  last  term  in  the  enetgy  equation  is 
a  triple  correlation  and  should  be  sisal  1  compared  to  double  correla¬ 
tions.  In  addition*  it  is  assumed  that  T  )^c(£ 

In  the  mean  steady  state,  then,  only  the  first  two  terms  in  bracket 
of  the  conservation  equations-  plus  the  first  terms  of  the  expansions 
are  retained  The  equations  are  now  prepared  for  introduction  of 
the  usual  boundary  layer  assumptions  which  lead  to  the  equations  of 
Section  2. 
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APPENDIX  B 

DISCUSSION  OF  THE  NUMBER  OF  INDEPENDENT  RELATIONS 
FOR  MOLECULAR  SPECIES 


Although  It  may  be  Intuitively  obvious  that  the  atomic 
mass  fractions, TV*.,  plus  the  rates  of  reaction,  UT{  ,  are  sufficient 
to  determine  all  of  the  luuleculat  mass  fractions,^;  ,  it  is  of 
interest  to  show  this  in  a  formal  way  and  also  to  show  how  many  of 
these  quantities  are  required.  Suppose  it  is  assumed  that  the 
system  contains  p  molecular  compounds  formed  by  reactions  between 
air  and  the  surface  material  and  that  these  molecular  species  are 
built  from  S  atomic  species.  Then  the  atomic  mass  fractions  are 
given  hv: 


where  is  the  number  of  atoms  of  atomic  species  •'***  contained 

in  molecular  species  L  end  may  be  zero.  From  the  theory  n:  linear 
equations  it  is  known  that  if  the  rank  of  the  matrix  of  «*,«.  is /L 
then  equation  (B-i)  represents  -A-  independent  equations  for  the  j5 
molecular  mass  fractions  rQ  .  In  order  for  equation  (B- 1 )  to  have 
any  solutions  at  all  it  is  required  that  there  be  ( S->t.  )  relations 
among  the  S  atomic  mass  fractions  TW  which  can  be^Obtained  by 
making  the  TWs  orthogonal  to  the  solutions,  if  any,  of  the  homo¬ 
geneous  equation  formed  with  the  transposed  matrix.  These  relations 
eliminate  any  superfluous  equations  from  the  syatcm.  In  addition,  it 
can  be  shown  that  the  races  of  reaction,  (A1^  ,  satisfy: 
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(B-2) 


equation  (B-2>  has  (  ? -rt* )  Independent:  scl.it  ions  so  that 
if  (a)  of  the  are  known  the  rest  can  he  determined.  These 

(  ?  )  production  rates  are  related  In  a  nnnl  Inrir  way  t"  tl*>  molecular 
mass  fractions  tnrough  tire  law  of  mass  action: 


Wc 


(B-3) 


The  products  arc  taken  over  all  species  which  enter  into  a 
given  reaction,  Ji  ,  while  the  sums  are  taken  over  all  of  the  L. 
reactions  which  occur.  The  quantities  and  are  the  forward 

and  backward  rate  coefficients,  respectively.  The  (p-.q,)  indepen¬ 
dent  equations  (B-3)  plus  the  /i.  Independent  equations  (B-l) 
represent  a  total  of  ?  equations  for  the  p  unknowns  rf;  .  If  Is 
assumed  that  among  the  roots  to  this  nonlinear  set  of  equation* 
there  is  one  real  root  for  each  between  zero  and  unity. 


In  addition  to  the  information  already  obtained,  it  is 
possible  to  determine  how  many  reactions  can  occur  between  th® 
chosen  set  of  molecules  and  to  find  a  set  of  reactions  which  can  be 
regrouped  to  express  the  actual  reactions.  The  molecular  weight  of 
molecular  species  t,  can  be  obtained  by  summing  the  weights  of  all 
atomic  species  which  it  contains: 


S 

£ 


(B-4) 
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Since  C<i^n  is  just  the  transposed  matrix  of  ,  it  is 

also  of  rank  n.  Therefore,  in  order  for  the  set  of  equations 
(B-A)  to  he  consistent  ft  is  necessary  that  the  satisfy  (f  - /i  ) 

relationships  obtained  in  a  scanner  similar  to  that  described  for  the 
T'«i'S  .  These  (P'-'t  )  relations  among  the  molecular  weights  can  be 
rearranged  to  give  the-  actual  reactions  themselves; 


•• 


a-*) 
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ESIIM'.I.’ON  OF  COMBUSTION  ZONE  THICKNESS 


An  estimate  of  the  thickness  of  the  combustion  zone  can  bt 
obtained  ill  the  following  iiuimer.  Let: 


], 


(C-I) 


U  -2> 


Differentiation  of  equation  (54)  with  respect  to  £  then  leads  to 
the  following  expression:  * 


f±itM ?  N  =  KV-f. 

v  Ai'Mf  )  I  ->-6* 
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From  the  del* It. 1 1 ion  of  J  hy  u^Uiuicm  (VJ)  the  change  in 

velocity  ratio  is  rclatrd  to  the  change  in  £  through: 

1^2  I  =  ^  (C-4) 


H*  £ j^v  and  are  the  same  order  of  twgnitude  it  seen:  reasonable 

t(»  define  a  characteristic  change  in  velocity  raliu  by  means  of: 


£>2*'-  ^  - 1 t- ^  1 

3  iis* 


(C-5) 


As  discussed  i"  Section  4  n't  i  =.  often  very  small  at  the 
combustion  zone.  Under  these  conditions  the  velocity  change  across 
the  combustion  zone  is  very  smalt.  In  the  tully  turbulent  region  of 
the  boundary  layer,  however,  the  velocity  changes  very  slowly  with 
distance  from  the  surface  so  that  in  the  turbulent  region  the  com¬ 
bustion  zone  nay  appear  rather  thick  even  when  the  velocity  change 
is  small .  in  the  sublayer  where  the  velocity  changes  extremely 
rapidly  the  zone  must  appear  thin.  When  *M  small  the 

boundary  layer  is  rather  sharply  divided  into  two  separate  regions. 
On  one  side  of  the  combustion  zone  oxygen  is  mixed  with  products  and 
inerts  while  on  the  other  side  gaseous  fuel  is  mixed  with  products 
and  inerts.  Since  depends  on  £  *,  temperature  and 

pressure  the  thickness  of  the  combustion  zone  car.  be  estimated  as  a 
function  of  B.  temperature,  and  pressure  at  the  comhusi.ion  zone  for 
any  particular  chemical  system. 
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FIG.  1  EFFECT  OF  FLOW  PAkAMETER,  ACCOMMODATION 
COEFFICIENT,  ,  ND  TEMPERATURE  ON  B 
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C  Qft,w  -  BX^  -  K«gJ  (BTU,' !  b) 
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FIG.  3  EFFECT  OF  FLOW  PARAMETER,  ACCOMMODATION  COEFFICIENT, 
AND  TEMPERATURE  O.l  ENTHALPY  POTENTIAL 
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WALL  TEMPERATURE  (°K) 


FIG.  5  EFFECT  OF  PRESSURE  AND  TEMPERATURE  ON  B 


TIC.  6  EFFECT  OF  PRESSURE  AMD  TEMPERATURE  OK 


FIG.  7  EFFECT  OF  PRESSURE  AND  TEMPERATURE 
OH  ENTHAlPY  POTENTIAL 
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FIG.  3  UNIVERSAL  ?LOT  FOR  INTEGRAND  OF 
(GRAPMITE  IN  AIR) 


FIG.  9  VARIATION  OF  WITH  MACH  NUMBER,  ALTITUDE 
AND  ANGULAR  LOCATION  ON  A  HEMISPHERE 
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PIG.  II  HFAT  FLUX  DISTRIBUTION  ON  HEMISPHERE- 
100,  00t  FT  ALTITUDE 
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ANUULAk  UiblANUE  lii'ir  STAGnATIGW  POIru ,  &  , 
(DEGREES) 


FIG.  13  DI  RIBUTIOH  OF  PAAAFETER  B  ON  HEMISPHERE-' 
100,003  FT  ALTITUDE 
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FIG.  14  DISTRIBUTION  OF  PARkMETER  0*.w  ON  HEM  SPHERE - 
100 ,0(X  FT  ALTITUDE 


-  66  - 


FIG.  15  HEAT  FLUX  DISTRIBUTION  ON  HEMISPHERE— 
50  ,00v.  FT  ALTITUDE 
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A-T.TITIT.E  ■  50,000  ft 


FIG.  16  MASS  FLUX  DISTR1 BUT ION  ON  HEMISPHERE-* 
50,000  FI  ALTITUDE 
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PIG.  19  DISTRIBUTION  OF  li»Ti5CML  ^  ON  HEMISPHERE 
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